Theory of Current and Shot Noise Spectroscopy in Single-Molecular Quantum Dots 

with Phonon Mode 
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Using the Keldysh nonequilibrium Green function technique, we study the current and shot noise 
spectroscopy of a single molecular quantum dot coupled to a local phonon mode. It is found that in 
the presence of electron-phonon coupling, in addition to the resonant peak associated with the single 
level of the dot, satellite peaks with the separation set by the frequency of phonon mode appear 
in the differential conductance. In the "single level" resonant tunneling region, the differential shot 
noise power exhibit two split peaks. However, only single peaks show up in the "phonon assisted" 
resonant-tunneling region. An experimental setup to test these predictions is also proposed. 

PACS numbers: 73.40.Gk, 72.70.+m, 73.63.Kv, 85.65.+h 
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Recent progress in engineering and fabrication of 
nanoscale electronic devices has made it possible to study 
the transport properties of molecular devices [|, the 
linear dimension of which is at least an order smaller than 
semiconductor quantum dots. In addition to their poten- 
tial industrial application, these devices provide an ideal 
test ground for the study of basic physics including the 
quantum size and many-body effects. The mechanism for 
electron conduction on such a small scale is not well un- 
derstood yet. However, the evidence for quantum nature 
of transport properties has been observed in differential 
conductance for molecular wires by several experimental 
groups § g, g, @, |, |, 0, 0, [if]. Theoretically, a lot of 
effort has been focused on the study of current-voltage 
characteristics of molecular wires based on either semi- 
empirical £jj, H or ab initio fl| ||, g Q 
methods. In contrast to semiconductor quantum dots, 
which is quite rigid in space, molecules involved in the 
electron tunneling process naturally possess the vibra- 
tional degrees of freedom which will inevitably react to 
the electron transfer through the molecules. So far, the 
influence of inelastic scattering process on the current- 
voltage characteristics of molecular wires has been con- 
sidered ]l7t while the influence on the current and its 
fluctuation (i.e., shot noise) has not been addressed for 
molecular dots. In this Letter, we use the Keldysh 
nonequilibrium Green function technique to calculate the 
current and shot noise through a single molecular quan- 
tum dot, for the first time. We focus on the effect of 
inelastic scattering process. A simplest Holstein-type 
model with a single phonon mode is employed to ad- 
dress the vibrational degrees of freedom in the molecular 
dot. All other complexity of real molecular devices, apart 
from interaction with phonon mode, is ignored. We find 
that in addition to the resonant peak associated with the 
single level of the dot, satellite peaks with the separation 
set by the frequency of phonon mode appear in the differ- 
ential conductance. In the "single level" resonant tunnel- 
ing region, the differential shot noise power exhibits two 
split peaks. However, only single peaks show up in the 




FIG. 1: Schematic diagram of a molecular quantum dot sys- 
tem. The dot is connected to two leads with couplings V L 
and T R . The dot electrons are also interacting with a single 
phonon mode. A gate electrode is capacitively attached to 
the dot to tune the single energy level. 



"phonon assisted" resonant-tunneling region. Therefore, 
due to the quantum nature of tunneling, the Fano factor 
is dramatically different from the Poisson limit even in 
the presence of inelastic process. These results may also 
be applied to the STM-bascd inelastic tunneling spec- 
troscopy around a local vibrational mode on surfaces p2[ . 

The model system under consideration is illustrated in 
Fig. [I]. It consists of a quantum dot connected with two 
normal conducting leads. The electrons on the dot is also 
coupled to a single phonon mode. The single level of the 
dot is tuned by a gate voltage. The system Hamiltonian 
is written as: 



H = Hr 
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The first two terms are respectively the Hamiltonian for 
electrons in the left and right non-interacting metallic 
leads: 
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where we have denoted the electron creation (annihila- 
tion) operators in the leads by c\ a (cka)- The quantity 
k is the momentum and a is the spin index and e k is the 
single particle energy of conduction electrons. The third 
term describes the phonon mode: 
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where ujq is the frequency of the single phonon mode, 
(a) is phonon creation (annihilation) operator. The 
fourth term describes the electron in the quantum dot: 

H D =J2[eo + Ka + ^)]dld a , (4) 

a 

where <v a (d a ) are the creation (annihilation) operators 
of dot electrons, e is the single energy level of the dot, 
A is the coupling constant between the dot electrons and 
the phonon mode. The last term represents the coupling 
of the dot to the leads: 

H T = [Vk*, a c{ a d a + H.c] , (5) 

where the tunneling matrix elements Vk a , a transfer elec- 
trons through an insulating barrier out of the dot. 

The current operators from the leads to the dot are, 
respectively, defined as: 

Il(r) = ie (Vka, a cl a d a - H.c.) . (6) 

keL(R);cr,a 

The terminal current is given by / = ((II) — (Ir))/2, 
where (...) denotes the statistical average of physical 
observables. Using the Keldysh nonequilibrium Green 
function formalism [ p3[ p4j , the terminal current can be 
obtained as |^5|, : 

I = de{Tr[(f L r L - f R r R )(G r - G a )} 

+Tr[(r L -r fl )G<]} , (7) 

where /l(r) are the Fermi distribution function of the left 
and right leads, which has different chemical potential 
upon a voltage bias /il — (J-r = eV, the coupling of the 
dot to the leads is characterized by the parameter: 

r L A R) = 2tt £ PH*),* ( £ ) ( 8 ) 

a 

with Pl(r),<j the spin-cr band density of states in the 
two leads, G^) and G^ a ,, are the Fourier transform 
of the dot electron retarded (advanced), G r J^)(t,t') = 
^fi6(±t^ft')({d a (t),dl t ,(t')}) i and lesser Green function, 

G aa' (*>*') = ^LCO^aC*))- Notc that the boldface no- 
tation indicates that the level-width function T and the 
dot electron Green functions are matrices in the spin 
space of dot electrons. 

The electrical current through a device is fluctuating 
in time even under a dc bias |^7j . Its fluctuation is char- 
acterized by the spectral density, which is given by the 
Fourier transform S(u>) of the current correlation func- 
tion g§: 

S(t,t') = ±({Ai L (t)ML(t')}), (9) 



where AIi,(t) = Ih{t) — {Ih{t))- The quantum statis- 
tical (non-equilibrium) average involved in the current 
correlation can be evaluated conveniently in the Keldysh 
Green function formalism. Following a tedious but stan- 
dard procedure, we obtain the spectral density of shot 
noise in the zero-frequency limit: 

S(u> - 0) = £ / de{-f L (l - /L)(Tr[(r L G' r ) 2 ] 
+Tr[(T L G a ) 2 ]) + z/ L Tr[r L G>] - i(l - / L )Tr[T L G<] 
+/ L Tr[r i G > r L (G r - G a )] + Tr[r i G > r L G < ] 

-(1 - / L )Tr[r L (G r - G a )r i G < ]} , (10) 

where G^ Q , is the Fourier transform of the greater Green 
function, G> a ,(t,t') = -i(d a (t)d f a ,(t')) . Equation © is 
a basic formula for the low-frequency shot noise through 
a quantum dot, which can take into account the many- 
body effects conveniently. It expresses the fluctuations 
of current through the quantum dot, an interacting re- 
gion, in terms of the distribution functions in the leads 
and local properties of the quantum dot, such as the oc- 
cupation and density of states. This formalism can be 
viewed as a generalized version of the two-terminal shot- 
noise formula for the non- interacting case [^9[ Q , which 
will become clear in the following discussion. It can also 
be used to study the current fluctuation in the case of 
spin-dependent transport, which is of much interest in 
the study of spintronics and single spin detection. 

Once the dot electron Green functions are known, the 
electrical current and shot noise can be calculated using 
Eqs. (g) and @. In the following, we calculate these 
Green functions, which should be carried out in the pres- 
ence of the leads. By performing a canonical transfor- 
mation, H = e s He~ s _with_5 = _(A/w ) Y, a Pa d <x{^ - 
a) one can obtain H = H e [+H p h- Here H e i = Hl + 
He + H D + H T and H ph = H x , where H L{R ^ = H L{R) , 

Hd = £J e o - X 2 /uj Q )dld a , H T « H t , and H x = H x . 
Consequently, the original dot-electron Green function 
can be decoupled as: 

G^W = Ti0(±t){{d a (t)Jl(0)}) el (X(t)XH0)) ph 
= G r J:)(t)(X(t)Xl(0)) ph , (11) 

where d a (t) = e iS ° lt d a e- iS ° lt , and X(t) = 
e iH P Ht Xe -iH ph t with x = exp [-(A/a;o)(a t - a)]. The 
renormalization factor due to the electron-phonon inter- 
action is evaluated to be f|: {X(t)X^(0)) ph = e~*W, 
where $(t) = (X/cj ) 2 [N ph (l - e^') + {N ph + 1)(1 - 
e-^"*)] with N ph = [exp(/3w ) - l]" 1 - We then apply 
the equation-of-motion approach to compute the Green 
function G^y . It consists of differentiating the Green 
function with respect to time, thereby generating 
higher-order Green functions which eventually have to 
close (by truncation in the presence of electron-electron 
interaction) the equation for the Green function G r ^ a \ 
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A little algebra gives rise to the Fourier transform: 



- (co - A) - S 



r(o) ' 



(12) 



where the energy shift A = \ 2 /u>o and the retarded (ad- 
vanced) self-energy due to the tunneling into the electri- 
cal leads are given by: 
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The full width of the resonance is then just the sum 
of elastic couplings to the two leads which is given by 
Eq. (|J). The Fourier transform of the full Green func- 
tion given by Eq. fxTu) can be obtained as: 



r> r ( a ) 



-g(2N ph + l) 



M2gjN ph (N ph + i)) 



i= 



w - (e - A) - Iujq - 



(14) 



where Ji(x) are the Bessel functions of complex argu- 
ment, the parameter g = (A/wo) 2 - The Green func- 
tion G^*, cannot be obtained by directly using the 
above equation-of-motion approach. By instead ap- 
plying the operational rules as given by Langreth Q 
to the Dyson equation for the contour-ordered Green 
function, one can show the following Keldysh equation 
for the lesser and greater functions |p6| : G < ^(cu) = 
G r (o;)Ey (Lu)G a (u>) . In the weak electron-phonon 
coupling limit as we are considering here, the contribu- 
tion to the self-energy from the electron-phonon 
interaction is negligible such that can be approxi- 
mated by the lesser (greater) self-energy due to the tun- 



neling into the two leads: X' 



i[f L T L + f R T R ] , and 



£> = -i[(l-/i)r- t + (l-/ii)r H ] . This relation enables 
us to rewrite the current and shot noise formulae as: 



Ji?jTr[T] 



(15) 



and 



s = ij de ^(i - h) + fnii - f R )mT] 

+ (/L-/«) 2 Tr[(l-T)T]}, (16) 

where we have defined the transmission coefficient ma- 
trix T = G a T L G r T R . We remark that for the nonin- 
teracting case (no electron-phonon and electron-electron 
interactions), the above two expressions are exact. These 
are just the Landauer-Buttiker formalisms developed for 
the noninteracting electron transport based on the scat- 
tering matrix theory. The connection between the two 
formalisms for the current was first established by Meir 
and Wingreen p3]. 




FIG. 2: Differential conductance, in units of 2e 2 /h, through 
the molecular quantum dot as a function of the Fermi energy 
Ep measured relative to the single level of the quantum dot 
in the presence (red-solid line, A = 0.4ljo) and absence (blue- 
dashed line, A = 0) of electron-phonon coupling. The energy 
is measured in units of the frequency of the phonon mode ujq , 
and r L = F R = 0.2o;o. 




FIG. 3: Differential shot noise power, in units of 2e 3 /h, 
through the molecular quantum dot as a function of the Fermi 
energy Ef measured relative to the single level of the quantum 
dot in the presence (red-solid line, A = QAujq) and absence 
(blue-dashed line, A = 0) of electron-phonon coupling. The 
energy is measured in units of the frequency of the phonon 



mode wo, and Y 



Q.2uj . 



For simplicity, we consider the coupling of the dot to 
the two leads to be symmetric T L = T R , and assume 
that the leads have broad and flat density of states (i.e., 
the wideband limit). For this case, the elastic couplings, 
from which the self-energy (|l^) can be determined, are 
independent of energy. In Figs. || and |^, we plot the 
zero-temperature differential conductance and shot noise 
power as a function of the Fermi energy (that is, the 
equilibrium chemical potential in the leads) measured rel- 
ative to the single level cq in the presence (with typical 
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value A = OAu>o, red-solid lines in the figures) and ab- 
sence (A = 0, blue-dashed lines) of electron-phonon cou- 
pling. Here all energy quantities are measured in units 
of the phonon mode frequency loq, and the elastic cou- 
plings r L = T R are fixed at 0.2wo- When the electron 
is not coupled to the phonon mode, only one resonant 
conductance peak shows up as the Fermi energy matches 
the single level in the quantum dot. In the presence of 
electron-phonon coupling, the overall spectrum is shifted 
by A = \ 2 /ojq. In addition to the main peak related 
to the single level, new satellite resonant peaks appear 
at the positive energy side. The separation between the 
conductance peaks is set by the frequency of the phonon 
mode ujq. Since at zero temperature no phonon modes 
are excited on the quantum dot, the electrons tunneling 
onto the dot can only excite phonon modes, which ex- 
plains why the satellite peaks are located at the positive 
energy region. It can be expected that there will appear 
satellite peaks at negative energy region at finite temper- 
atures, where thermally excited phonons can be absorbed 
by the electrons tunneling onto the dot. Moreover, the 
intensity of the satellite peaks is much smaller than the 
main resonant peak because they are evolving from the 
emission of phonon modes, which is controlled by the 
electron-phonon coupling. Correspondingly, the differ- 
ential shot noise power exhibits two peaks located sym- 
metrically around the position where the conductance 
peak associated with the single level of the dot is located. 
The origin of this behavior is that no noise is generated 
when the transmission probability T = or 1, while it 
is generated maximally in between. In this "single level" 
resonant-tunneling region, the Fano factor defined as the 
ratio of shot noise to current is significantly enhanced. 
Since the transmission probability corresponding to those 
satellite conductance peaks (from the phonon emission 
process) is small, the differential noise power is approx- 
imately proportional to the conductance and only show 
a single peak rather than two peaks. Therefore, in the 
"phonon assisted" resonant-tunneling region, the Fano 
factor is smaller than one, in contrast to the case of "sin- 
gle level" resonant tunneling. Therefore, due to the quan- 
tum tunneling nature, the Fano factor for the present sys- 
tem differs the Poisson limit (equal to 2), where electrons 
diffuse in an uncorrelated way. 

To experimentally investigate the current and shot 
noise spectroscopy of molecular quantum dots, we sug- 
gest a setup where a gate electrode is attached to the dot 
capacitivcly, as shown in Fig. |l|. A gate voltage is ap- 
plied solely to tune the level position on the dot while a 
bias voltage is applied in the linear response regime. The 
merit of the proposed setup is that the internal electronic 
and ionic charge buildup on the dot, which is generated 
by a nonlinear response, could be avoided safely. Exper- 
imental effort along this direction is in progress |33j . 

To summarize, using the Keldysh noncquilibrium 
Green function technique, we have studied in this work, 



to the best of our knowledge, for the first time the cur- 
rent and shot noise spectroscopy of a single molecular 
quantum dot coupled to a local phonon mode. We show 
that, in the presence of electron-phonon coupling, in ad- 
dition to the resonant peak associated with the single 
level of the dot, satellite peaks with the separation set 
by the frequency of phonon mode appear in the differen- 
tial conductance. In the "single level" resonant tunnel- 
ing region, the differential shot noise power exhibits two 
split peaks. However, only single peaks show up in the 
"phonon assisted" resonant-tunneling region. The Fano 
factor is different from the Poisson limit even in the pres- 
ence of inelastic process. The derived current and noise 
formalism can also be applied to more complicated sys- 
tems. 
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